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For humans, wheat is the most important source of calories, but it is also a source of 
antioxidant compounds that are involved in the prevention of chronic disease. Among the 
antioxidant compounds, phenolic acids have great potential to improve human health. 
In this paper we evaluate the effect of environmental and genetic factors on the 
phenolics content in the grain of a collection of tritordeums with different cytoplasm and 
chromosome substitutions. To this purpose, tritordeum flour was used for extraction of 
the free, conjugates and bound phenolic compounds. These phenolic compounds were 
identified and quantified by RP-HPLC and the results were analyzed by univariate and 
multivariate methods. This is the first study that describes the composition of phenolic 
acids of the amphiploid tritordeum. As in wheat, the predominant phenolic compound is 
ferulic acid. In tritordeum there is great variability for the content of phenolic compounds 
and the main factor which determines its content is the genotype followed by the 
environment, in this case included in the year factor. Phenolic acid content is associated 
with the substitution of chromosome DS1D(1H*) and DS2D(2H*), and the translocation 
1 RS/1 BL in tritordeum. The results show that there is high potential for further improving 
the quality and quantity of phenolics in tritordeum because this amphiploid shows high 
variability for the content of phenolic compounds. 
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1. INTRODUCTION 

Cereals, particularly wheat, are the most important food source 
for human. Numerous studies have demonstrated the key role of 
whole grains in human health benefits, since they are known to 
be protective against chronic diseases. These data provide further 
support for recommendations to increase consumption of whole 
grains and its use as healthy products (Marquardt et al., 2004; 
Flight and Clifton, 2006; de Munter et al, 2007; Jonnalagadda 
et al., 2011). A wide range of phytochemicals with antioxidant 
activity present in cereal grains has been indicated to be responsi- 
ble for these health protection qualities. These compounds belong 
to chemical groups such as polyphenols, carotenoids, and plant 
sterols. Among the different antioxidants present in wheat, phe- 
nolic compounds seem to have the greatest potential of being 
beneficial to health (Beta et al, 2005; Mpofu et al, 2006). 

Phenolic acids are the most frequent phenolic compounds in 
cereals; these can be found as free, bound and conjugated forms. 
Most phenolic acids are bound by ester-linked to the cell wall 
polymers (Irakli et al., 2012). The main phenolic acids in wheat 
are ferulic and p-coumaric acids, both associated with cell-wall 
constituents (Zhou et al, 2004, 2005; Okarter et al, 2010). It has 
been described that factors such as the genotype, environmen- 
tal factors and genotype-environment interactions can influence 
total phenolic compounds content and the relative content of 
different phenolics fractions. Yu and Zhou (2005) described the 
effect of the location, on the content of phenolic compounds; 
Moore et al. (2006) observed significant correlations between 



temperature stress or solar radiation and some antioxidant con- 
tents. Furthermore, Fernandez-Orozco et al. (2010) showed that 
environmental factors did not equally affect to all phenolic com- 
pounds fractions, the fractions of free phenolic compounds and 
conjugates being more affected than the fraction of bound phe- 
nolics. 

In addition, the content of phenolic compounds is affected by 
the genotype and the genotype x environment interaction. The 
results of Mpofu et al. (2006) showed that environmental effects 
on the content of phenolic compounds were considerably larger 
than the genotype effect, and genotype x environment inter- 
action was small for all parameters as compared to both main 
factors . The study of Fernandez-Orozco et al. (2010) also showed 
that the environmental effects were larger than genotypic differ- 
ences, especially for the low abundance of free and conjugated 
phenolic acid fractions. Li et al. (2008) reported differences in 
phenolic content between varieties, and they suggested that it may 
be possible to develop varieties with a high content of phenolic 
acids. The results of Shewry et al. (2010) corroborated previous 
reports and demonstrated that heritable variation in the con- 
tent of bioactive components can be exploited by breeders to 
develop new cultivars with enhanced health benefits. Triticeae 
wild species might also be a valuable source of variability, but 
transferring this variability to cultivated species presents multi- 
ple problems. That is not the case of the wild species Hordeum 
chilense Roem. et Schultz. whose variability for traits of inter- 
est may be used in both durum and bread wheat breeding using 
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the amphiploids called tritordeums ( x Tritordeum Ascherson 
et Graebner) (Martin and Sanchez-Mongelaguna, 1982; Martin 
et al, 1998, 1999). 

Moreover, other qualities of tritordeum confer on it the 
potential to be part of a healthy diet and it has recently been 
introduced to the market as an alternative food cereal (www. 
tritordeum.com). For example, tritordeum has a high con- 
tent of carotenoids, such as lutein, which has been associated 
with prevention against different diseases (Atienza et al., 2007; 
Mellado-Ortega and Hornero-Mendez, 2012). To the present, 
little is known about the content and variability of phenolic 
compounds, especially phenolic acids of tritordeum. In addi- 
tion, factors affecting the content of phenolic compounds have 
not been determined and it is not known whether they are 
the same as those involved in the content of phenolic com- 
pounds in wheat. Furthermore, little is known about the genetic 
and chromosomal location of the trait of high/low content of 
phenolic compounds. In order to carry out an efficient improve- 
ment of this character to increase health benefits, more infor- 
mation on the factors affecting phenolic content needs to be 
collected. 

The aim of the present work was to evaluate the whole grain 
phenolic acids content of accessions of hexaploid tritordeum and 
to assess the effect of the environment and some genetic factors 
such as the presence of barley or wheat cytoplasm; substitu- 
tion of chromosomes IH'^'', 2H^'', and 5H"'^; and the IRS/IBL 
translocation, on phenolic acid content. 

2. MATERIALS AND METHODS 

2.1. PLANT MATERIAL 

Thirty seven accessions of hexaploid tritordeum originating 
from different stages of the tritordeum breeding program 
were evaluated (Table 1). The previous genotyping of the tri- 
tordeum lines (Castillo et al., 2013) allowed their characteriza- 
tion. Different groups were identified according their genomic 
constitution (1) cytoplasm, (2) chromosome substitutions, and 
(3) IRS/lBL translocation. Tritordeum grains were ground using 
a cyclonic mill (Ciclotec lab mill, Tecator) with 1 mm mesh 
to produce wholemeal flour. Samples were stored at 4°C until 
analysis. 

2.2. PHENOLIC COMPOUNDS EXTRACTION 
2.2. 1. Extraction of free ptienolics 

Free phenolics (and also bound and conjugate phenolics) were 
extracted using the method previously reported by Okarter et al. 
(2010) with some modifications. One hundred milligram of 
wholemeal flour was blended with 670 |il of 80% chilled ace- 
tone in 2 mL tube. Inside of the tube we added a 5 mm stainless 
steel ball to facilitate the mixing in a Retsch MM200 ball mill for 
1 min at 60% of power. The mixture was centrifiiged at 2500 x g 
for 10 min The supernatant was removed and the remaining 
pellet was extracted two times more with 80% chilled acetone 
following the same process. The supernatants were pooled and 
evaporated in a Speed Vac vacuum centrifuge (Savant) to dryness. 
The solution was then reconstituted with methanol/hydrochloric 
acid (1 M; 85:15, v/v), filtered through a 0.45 |xm filter, and stored 
at — 20°C until analysis. 



Table 1 | Tritordeum collection used in this study. 



Line 


Year 


Cytoplasm 


IRS/1 BL 


Substitutions 


223 


10 


11 


12 


n. cnilense 


N 


N 


294 


10 


11 


12 


Wheat 


N 


C -1 / 1 LI ch \ 


295 


10 


11 


12 


n. cnilense 


N 


N 


296 


10 


11 


12 


Wheat 


N 


N 


320 


10 


11 


12 


n. cnilense 


N 


N 


322 


10 


11 


12 


n. cnilense 


N 


N 


323 


10 


11 


12 


n. cnilense 


N 


N 


324 


10 


11 


12 


n. cnilense 


N 


N 


325 


10 


11 


12 


n. cnilense 


N 


N 


326 


10 


11 


12 


n. cnilense 


N 


N 


327 


10 


11 


12 


n. cnilense 


N 


N 


328 


10 


11 


12 


Wheat 


N 


c 1 Pi / 1 LI ch \ 


400 


10 


11 


12 


Wheat 


N 


DS2D(2H ) 


409 


10 


11 


12 


vvneat 


D 

n 


Pi C C Pi / C LI ch \ 


410 


10 


11 


12 


H. chilense 


D 

n 


M 

N 


411 


10 


11 


12 


H. chilense 


D 

n 


M 


412 


10 


11 


12 


H. chilense 


N 


N 


413 


10 


11 


12 


H. chilense 


R 


N 


414 


10 


11 


12 


vvneat 


D 

n 


M 

In 


415 


10 


11 


12 


Wheat 


R 


ps c o pi / o LI ch \ 


416 


10 


11 


12 


Wheat 


R 


Pi CI Pi /o Lich \ 


417 


10 


11 


12 


H. chilense 


D 

n 


UozUlzn ) 


418 


10 


11 


12 


M A 


M 


M 


419 


10 


11 


12 


Wheat 


N 


N 


421 


10 


11 


12 


Wheat 


R 


r\cor\/oi_iGh\.p\ccp\/ci_ichn 


422 


10 


11 


12 


Wheat 


N 


N 


423 


10 


11 


12 


H. chilense 


N 


N 


424 


10 


11 


12 


Wheat 


R 


Db2D(2H );DS5D(5H ) 


425 


10 


11 


12 


Wheat 


R 


UbzU(zH ),IJbblJ(bn ) 


427 


10 


11 


12 


Wheat 


N 


N 


429 


10 


11 


12 


Wheat 


R 


N 


430 


10 


11 


12 


H. chilense 


R 


DS5D(5H*) 


431 


10 


11 


12 


H. chilense 


N 


DS5D(5H*) 


432 


10 


11 


12 


H. chilense 


R 


N 


436 


10 


11 


12 


H. chilense 


N 


DS5D(5H*) 


621 


10 


11 


12 


NA 


NA 


NA 


631 


10 


11 


12 


Wheat 


N 


N 



10:11:12, years 2010, 2011 and 2012 respectively: W, no translocation or no 
substitutions: Fl, 1RS/1BL translocation: NA no available. 

2.2.2. Extraction of bound pfienolics 

Insoluble-bound phenolics were extracted from the residue after 
the free phenolics extraction. The residue was then washed three 
times with distilled water to remove the water soluble sugars 
which may gelate and hinder the extraction of bound phenolics. 
After washing, the residual pellet was digested with 1.5 ml of a 2N 
sodium hydroxide solution. The mixture was blended in a Retsch 
MM200 ball mill for 1 min at 60% of power and incubated for 
1 h at room temperature under gentle shaking in darkness. Both 
the sodium hydroxide and the mixture of sodium hydroxide and 
the pellet were degassed with nitrogen. After digestion, the tubes 
were centrifuged for 2 min at 8000 x g and each supernatant was 
transferred to a new tube. The digestion was stopped with fumant 
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FIGURE 1 I RP-HPLC chromatograms of the bound (A), conjugated (B), and free (C) phenolics fractions from tritordeum grain. 



chloride acid (75%) to reach a pH 2. We extracted the pheno- 
hcs by adding 670 |il of ethyl acetate, mixing by vortex for 1 min 
and centrifuging for 10 min at 7000 x g. The supernatant was 
transferred to a clean tube. This process was repeated three times 
and the supernatants pooled. The ethyl acetate was evaporated 
in a Speed Vac vacuum centrifuge (Savant). The solution was 
then reconstituted with mehtanol/hydrochloric acid (1 M; 85:15, 
v/v), filtered through a 0.45 \im filter, and stored at — 20°C until 
analysis. 

2.2.3. Extraction of soluble-conjugated phenolic compounds 

Soluble-conjugated phenolics were extracted from half volume of 
free phenolics extracts using the method reported above for the 
bound phenolics with some minor modifications. To the residue 
from evaporation of the acetone of half volume of free phenolics 
extracts, we added 1.5 ml of 2N sodium hydroxide and from this 
step we followed the same process like in the extraction of bound 
phenolics. 

2.3. HPLC METHOD 

The phenolics extracts (100 |il) were applied to a Eclipse 
XDB-C18 reverse phase analytical column (4.6 x 150 mm, 



5 [xm particle size; Agilent Technologies) using a 1200 
Series Quaternary LC System liquid chromatograph (Agilent 
Technologies) with a DAD UV-V detector. We eluted the phe- 
nolic compounds using as solvent A a mix of water and 0.01% 
trifluoroacetic (TFA), and as solvent B a mix of acetonitrile 
(ACN) and 0.01% of TFA, according to the following method: 
1ml min^', and a gradient from 7% to 16% of solvent B (to 
100% with solvent A) over 22.5 min followed by a 65% of solvent 
B over 5 min to wash the column. Absorbance was monitored 
with the DAD UV-V module at 280 nm. 

2.4. EXPERIMENTAL DESIGN AND STATISTICAL ANALYSIS 

All analyses were conducted with the statistical software R ver- 
sion 2.12.1 (R Development Core Team, 2012). The experimental 
design was a completely randomized design replicated for 3 years 
(2010-2012) with a 2m^ plots. Data were adjusted to a lineal 
model with the function Im and factors effects were checked by 
a analysis of variance with the function anova. The lineal mod- 
els were adjusted with the factors "Year," "Line," "Substitution," 
"IRS/IBL" and "Cytoplasm" (see models formula in section 3). 
Factor "Year" has three levels (years 2010, 2011, and 2012), 
factor "Line" has 37 (see Table 3), factor "Substitution" has five 
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Table 2 | Table with the mean percentage and standard deviation of 
each phenolic acid relative to the total phenolic compounds. 



Phenolic compound 


Mean (%) 


SD 


Bound fsrulic acid 


ocj.yu 


K c^n 


Bound vanillic acid 


u.uu 


n nn 
u.uu 


Bound p-coumaric acid 


J.DU 


1 in 

1. lU 


Bound syringic acid 


U. lU 


n in 


Bound caffeic acid 


0.10 


0.10 


Bound dimsrs 


o. /u 


n on 

U.oU 


Fr66 fcrulic acid 


1 .oU 


1 on 

I.UU 


FrG6 vanillic acid 


u.ou 


n AO 

U.^U 


FrsG p-counnaric acid 


n 90 

U.oU 


1 in 

1. 1 U 


FrsG syrinQic acid 


1 R Rn 

1 O.DU 


R QO 

u.yu 


FrsG caffsic acid 


n AO 


1 .DU 


FrGG dimGTS 


/.DU 


1 p^n 

I.OU 


Conjugatod fGrulic acid 


D.ZU 


o in 

z. lU 


Conjugatod vanillic acid 


u. / u 


n on 
u.zu 


Conjugatod p-coumaric acid 


n 7n 


n Rn 

U.DU 


Conjugatod syringic acid 


n Qn 

u.yu 


n Qn 
U.oU 


Conjugatod caffoic acid 


n Kn 


n Qn 
U.oU 


Conjugatod dimGrs 


1 Qn 
1 .yu 


n Qn 

U.oU 


lULdl IcIUIlL. dL-IU 


76 40 


6 40 


Total vanillic acid 


1.50 


0.60 


Total p-coumaric acid 


4.50 


1.70 


Total syringic acid 


16.60 


6.90 


Total caffoic acid 


1.00 


1.60 


Total bound 


72.60 


6.60 


Total froG 


18.50 


6.60 


Total conjugatGd 


8.90 


2.90 


Total dimGrs 


16.20 


2.00 



7b calculate the mean and standard deviation values shown in the table we used 
all the lines and years, sd, standard deviation; dimers, dimers of phenolic com- 
pound; "total" is the sum of bound, free and conjugated phenoiics. Percentages 
are calculate using absorbance data. 
SD: standard deviation en %. 



levels (DS1D(1H^'^), DS2D(2ff'^), DS5D(5H*), DS2D(2ff'^) 
DS5D(5H'^'^) and N), "IRS/IBL" has two levels (R (IRS/BL 
translocation) and N (no translocation)) and "Cytoplasm" has 
two levels (cytoplasm from H. chilense and from Wheat). 

The normality and heteroscedasticity assumptions were tested 
by plotting the residuals versus the predicted values and Q-Q 
plots. In the cases where the conditions of data normality and 
homogeneity of variances were violated the Box-Cox transfor- 
mation was applied (function powerTransform, package car) (Fox 
and Weisberg, 2011). The differences between tritordeum lines 
were assessed using post hoc multiple-comparison test (function 
glht, package multcomp) (Hothorn et al., 2008). The method 
used for multiple-comparison was a Tukey contrast with a 
p-values adjustment type free (see package multcomp for more 
information). 

To identify the relative importance of each factor in the 
total set of data, we performed a multivariate analysis of vari- 
ance (MANOVA) with the function Im. The multivariate models 
were adjusted with the factors "Year," "Line," "Substitution," 



"IRS/lBL" and "Cytoplasm" (see models formula in section 3). 
As an exploratory graphic, we generated a Principal Components 
Analysis (PCA) to demonstrate trends of tritordeum lines in 
relation to the variables studied using the function PCA (pack- 
age FactoMineR) (Husson et al., 2013). The package maptools 
was used to avoid the label overlapping in the ordination plot 
(Lewin-Koh et al., 2012). 

3. RESULTS 

3.1. IDENTIFICATION OF PHENOLIC COMPOUNDS 

For identification of phenolic compounds by RP-HPLC we 
used pure compounds from SIGMA-ALDRICH (http://www. 
sigmaaldrich.com/). The pure compounds ferulic acid (FA, cat# 
46278), p-coumaric acid (pCA, cat* C9008), caffeic acid(CA, 
cat# C0625), vanillic acid (VA, cat# 94770) and syringic acid 
(SA, cat# S6881) were dissolved in methanol at a concentra- 
tion of 1 mg/ml and injected in HPLC according to the pro- 
gram described in section 2. We used the chromatograms of 
pure compounds for their identification in the grain samples 
analyzed (Figure SI). For the identification, we used mainly 
the retention times (ferulic acid ~17min; p-coumaric acid 
~15.2min; caffeic acid ~11.2min; vanillic acid ~10.4min.; 
syringic acid ~11.7min). In Figure SIA we showed the chro- 
matogram of a mix of the same amount of each phenolic 
compound. The area and height at 280 nm of the monomeric 
peaks were very similar. This means that the response factor 
(absorbance by weight unit) was very similar for all pure pheno- 
lic compounds studied in this paper. The peaks appearing from 
minute 19 onwards were considered dimeric phenoiics (phenolic 
compounds all together). We tried to identify the phenolic com- 
pounds peaks in the grain samples using retention time, but we 
found small variations in the compounds peak migrations times 
and also small unknown peaks near phenolic compounds peaks. 
In these cases, we used the spectrum of the pure compounds from 
210 to 400 nm and we compare it with the peaks present in the 
samples (Figure S2) using the tool for detection of pure spectra 
of the ChemStation software. When the peak spectrum showed 
no coincidence above 80% with any of the pure spectra no peak 
was identified. Moreover, we did not associate the peaks of the 
dimeric region with any particular phenolic compound and we 
could not measure the total amount of each individual pheno- 
lic compound. We extracted free, bound and soluble-conjugated 
phenolic compounds from whole flour of a set of Tritordeum lines 
assayed during three consecutive years. The different fractions of 
phenolic compounds were injected into liquid chromatograph 
columns as described in section 2 (Figure 1). In Figure lA we 
plot a typical chromatogram of the bound phenolic fraction. This 
chromatogram presented pCA, FA and dimers. In the fraction 
of soluble-conjugated phenoiics (Figure IB) we found all phe- 
nolic compounds studied in this work (FA, pCA, CA, VA, and 
SA), although some of them at very low amounts. Finally, the 
fraction of free phenoiics VA, SA, FA, and dimers is presented 
(Figure IC). 

3.2. QUANTIFICATION OF PHENOLIC ACID COMPOSITION 

The relative quantification is shown in Table 2. The values shown 
are the mean percentages and standard deviation of all lines and 
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FIGURE 2 I Principal components analysis. (A) Tlie quantitative 
variables are represented. (B) The individuals lines and the centroids of 
levels factors "Year," "Sustitutions," "1RS/RL," and "Cytoplasm" are 
represented. X "B.FA," "F.FA," "CPA," "Total. PA" are bound, free, 
conjugated and total ferulic acid; "B.VA," "PVA," "C.VA," "Total. VA" are 
bound, free, conjugated and total vanillic acid; "B.pCA," "PpCA," 



"C.pCA," "Total. pCA" are bound, free, conjugated and total p-coumaric 
acid; "B.CA," "PCA," "CCA," "Total. CA" are bound, free, conjugated 
and total caffeic acid; "B.SA," "PSA," "C.SA," "Total. SA" are bound, 
free, conjugated and total syringic acid; "Total. B," total bound phenolic; 
"Total. p" total free phenolic; "Total. C," total conjugated phenolic; 
"Total. phenolic," total phenolic; "Total. Dim," total dimer phenolic. 



years. The most abundant compound in the fraction of bound 
phenolics is ferulic acid followed by p-coumaric. Ferulic acid is 
also the main compound in the fraction of conjugated phenolics. 
In the fraction of free phenolics is more abundant syringic acid. 
In general, the main compound quantitatively is the ferulic acid 
and the second most abundant is the syringic acid. With respect 
to the fractions, the highest content of phenolic compounds are 
in the bound fraction and the lowest in the conjugated. 

3.2. 1. Exploratory multivariate analysis (MANOVA and PC) of 
phenolic compounds content 

We carried out a multivariate analyze to explore the whole data 
set. We fitted two models because the factor "Line" did not 
have repetitions inside of each year. Therefore, the Model 1 was: 
"variable ~ Line + Year". The analysis of variance table showed 
that "Line" and "Year" had a highly significant effect on the 
data set (Table SI; Model 1). The Model 2 was fitted with main 
factors Year, Substitutions, translocation IRS/IBL, Cytoplasm 
and the two interactions between them ('variable ~ Year + 
IRS/lBL + Substitutions-I- Cytoplasm + Year x IRS/lBL + 
Year x Substitutions + Year x Cytoplasm + IRS/lBL x 
Substitutions + IRS/lBL x Cytoplasm + Substitutions x 
Cytoplasm'). The main factors Year, IRS/lBL and Substitutions 
had a significant effect on the data set (Table 81; Model 2). The 
interactions Year x Substitutions and IRS/lBL x Substitutions 
also showed a significant effect. These results were used to sim- 
plify the models for the univariate analysis of each variable. 

Principal component analysis (PC), based on correlation 
matrix, was performed using the different phenolic compounds 
contents and the factors with significant effects. According to 



the Kaiser-Harris criterion (Kaiser, 1974) we selected eight com- 
ponents with eigenvalues greater than 1. We used the principal 
component analysis as an exploratory analysis and not as a tool 
to reduce the number of variables. For this reason, we plotted 
the variables and factors regarding the principal component one 
and two, which accumulate the higher variability. Most of the 
variables are correlated with principal component 1 (PCI) being 
positive for aU of them (Figure 2A). The PC2 showed a high 
correlation with the amounts of free vallinic acid, total vanil- 
lic acid, total free phenolic, free and total syringic acid. In this 
case, the amount of free and total syringic acid had a negative 
correlation with free vallinic acid, total vanillic acid and total 
free phenolics. The amount of free p-coumaric acid, free caf- 
feic acid, bound syringic acid, free ferulic acid, bound vallinic 
acid and bound syringic acid showed a low correlation with PCI 
and PC2. 

The factors are represented on the PCI and PC2 (Figure 2B). 
The factor "Year" had all its levels spread in the four differ- 
ent quadrants. The year 2011 and 2012 showed some variation 
along PC2; the year 2010 showed differences with respect to 
the year 2012 along of the PCI and with the 2011 in both 
(PCI and PC2). The "Substitutions" factor presented a wide 
dispersion of their levels. Disomic substitution DS2D(2H'^'') 
and DS5D(5H"^*^) showed the higher differences regarding to 
the "no substitution" level, mainly along PCI. Disomic substi- 
tution DS2D(2H'^'') showed the greatest differences regarding 
to the level "no substitutions" along both PCI and PC2. The 
two levels of the factor "IRS/lBL" translocation showed dif- 
ference in PCI. The "Line" factor showed a wide dispersion, 
although there is some accumulation in quadrant 3 (— , — ), where 
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Table 3 | Table with the percentage of variance explained by each 
factors, derived from the ANOVA table of Model 1 ('variable ~ Line + 
Year'). 



Model 1 




SS percentage 




Line 
(Df = 36) 


Year 
(Df = 2) 


Residuals 
(Df = 72) 


IOTBI TGrUllC 3CIQ 


oo.u 1 


O / .Z 1 


9Q 9fi 


ConjuQBtsd fsrulic scid 


9Q 7n 


9R 7"^* 
ZD. /o 


AA R7 


Fr66 fGrulic 3cid 


OD.Zlj 


n OA 


71 
/ 1 


RniinH foriilir' ar^iH 


36 02** 


32 36* 


31 63 


Tot3l vsnillic scid 


oU. / o 


9R Al* 


lA 


Conjug3tsd vanillic acid 


AO n7* 
^z.u / 


7 71 ** 
/ . / 1 


99 

uU. ZZ 


1 U Ldl p LUU 1 1 Idi lU dUlU 


31 11 * 


35 29* 


33 60 


•^Ul IJ Uyd LcU p CUUI 1 Id 1 IL. dL-IU 


40 79* 


QQ OA* 
oO.O'+ 
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Total free 
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Total conjugated 


34.06. 
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Total phenolic acids 
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Conjugated dimers 
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SS, sum of the squared; Df, degree of freedom. 
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the control levels of factors "Substitutions" and "IRS/IBL" are 
found. 

3.2.2. Univariate analysis ofplienolic compounds content 

We carried out a univariate analysis of variance to check which 
variables were affected by the factors studied. We fitted two 
models for the same reason that in the multivariate analysis of 
variance. Moreover, the MANOVA Model 1 and 2 showed factors 
with significant effects, consequently, we conducted univariate 
analysis of each of the variables using both models. In MANOVA 
Model 1 all factors, "Line" and "Year," were significant, for that 
reason we use the same model ('variable ~Line + Year'). In 
this analysis the factor "Line" had a significant effect (with a p- 
value <0.05) for all variables except conjugated FA, total VA, 
bound pCA, total CA and total conjugated (Table 3). The per- 
centage of variance explained by the factor "Line" range between 
57.6% and 28.2%. The factor "Year" had a significant effect on all 
variables except free FA, conjugated SA and free dimers. The per- 
centage of variance explained by the factor "Year" range between 
41.2% and 7.71%. 

From the results of MANOVA Model 2, we simplified the 
model by eliminating terms that were not significant. Therefore, 



we used the following Model 2 in the univariate analysis: "vari- 
able ~Year + Substitutions + IRS/lBL -|- Year x Substitutions 
-|- Substitutions x IRS/BL". In this model the factor "Year" 
showed behavior almost identical to that shown in Model 1, 
showing even a very similar explained variance (Table 4). The 
factor "Substitution" showed a significant effect for all variables 
except bound pCA, total SA, free SA, total bound and total free. 
The variance explained by the factor "Substitution" ranged from 
5.8% to 32.6% of the total variance. The factor "IRS/lBL" had a 
significant effect on total FA, conjugated FA, bound FA, bound 
pCA and total bound. The percentage of variance explained 
by this factor ranged from 2.2% to 3.6%. Interaction "Year x 
Substitutions" had significant effects on the variables conjugated 
FA, total conjugated and conjugated dimers, with percentages 
of explained variances from 15.8% to 18.7%. The last terms of 
Model 2, the interaction of "Substitution x IRS/lBL" showed 
a significant effect for the variables total FA, bound FA, conju- 
gated VA, total pCA, conjugated pCA, free pCA, conjugated SA 
and conjugated CA. The percentages of variance explained by this 
interaction are between 3.3% and 10.0%. The factor "IRS/lBL" as 
shown in Figure 3 increased the FA content in general,(soluble- 
conjugate, bound and total ferulic), bound pCA and total bound 
phenolics. 

For the factor "Substitution" we carried out a comparison 
of means between the different levels (Figure 4). Chromosome 
substitution DS1D(1H'^'') caused an increase of total pCA. The 
substitution DS2D(2H^'') was the substitution with the higher 
effect on phenolic compounds. This substitution increased total 
pCA, conjugated pCA, conjugated CA, total dimers and free 
dimers. In fact, three of the four lines with the highest contents 
of phenolic compounds had the substitution DS2D(2H'^'')(lines 
400,415, and 416). 

4. DISCUSSION 

Phenolic compounds of whole grain are an important trait 
because they have several human health benefits. Wild cereal 
species are an important source of variability that can be used 
by plant breeders to produce healthier wheat varieties. The 
amphiploid tritrodeum is an excellent way to access the vari- 
ability present in the wild species H. chilense. In this paper we 
characterize the phenolic acid composition of a collection of tri- 
tordeums, and we assess the effect of growing season, the barley or 
wheat cytoplasm, chromosomal substitutions and the widespread 
chromosomal translocation IRS/lBL. 

Until now, there was little information on the composition of 
phenolic acids, the main compounds with antioxidant activity, 
of the amphiploid tritordeum. Ferulic acid is the predominant 
phenolic acid found in whole grains (Krygier et al., 1982; Mattila 
et al, 2005) and whole wheat (Moore et al, 2005; Li et al, 2008). 
These authors have reported ferulic acid percentages to the total 
of phenolic compounds in the range of 72-85%, in the case of Li 
et al. (2008) and 83.5-89.5% in the case of Moore et al. (2005). 
The percentages of ferulic acid with respect to total phenolics 
in our tritordeum collection (on average 76.38%) are similar to 
those found in wheat. In our study, 90.2% of total ferulic acid is 
found in the fraction of bound phenolic acids. These values agree 
with those reported in wheat, since several authors have shown 
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Table 4 | Table with the percentage of variance explained by each factors, derived from the ANOVA table of Model 2 ('variable ~ Year + 
Substitutions + 1RS/1BL + Year x Substitutions + Substitutions x 1RS/BL'). 
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FIGURE 3 I Content of phenolic compounds with a significant effect of the translocation 1RS/RL (p-value < 0.05; Table 4). N, no translocation; R, 
translocation 1RS/RL. 
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FIGURE 4 I Phenolic add content and multiple-comparison results of 
the different substitution lines. Gray bars are significantly different 
(p-value < 0.05 ) compared to the no substitution (N). Comparison of 
means (Tukey contrast) between levels of factor "Substitutions" took 



place only when the ANOVA table of Model 2 (Table 4) showed a 
significant effect of this factor. Figure shows only those compounds 
where any substitution level showed differences with respect to no 
substitution level (N). 



percentages of ferulic acid in the bound fraction from 89% to 97% 
(Adorn et al., 2003; Moore et al., 2005; Li et al, 2008). 

Other phenolic acids were also found in significant quanti- 
ties in whole wheat. Okarter et al. (2010) have reported that the 
second most abundant phenolic acid found in whole wheat is p- 
coumaric acid. In this study, the second more abundant phenolic 
compound is the syringic acid and not the p-coumaric acid, which 
was found in smaller amounts than syringic acid (16.6% of total 
phenolic acids is syringic acid and only the 4.52% of total pheno- 
lics is p-coumaric acid). The higher ratio of syringic acid in grain 
has been also reported by Li et al. (2008) in spring and spelt wheat 
and some varieties of soft wheat (Moore et al., 2005). 

Syringic and vanillic acids are hydroxybenzoic acids, and fer- 
ulic acid, p-coumaric acid and caffeic are hydroxyninnamic. 
The hydroxyninnamic and hydroxybenzoic acids belong to dif- 
ferent biosynthetic pathways; hydroxybenzoic acids belong to 
aminobenzoate pathway, and the hydroxyninnamic acids to the 
phenylpropanoid biosynthesis pathway. Figure 2A shows a high 



correlation of the total hydroxybenzoic acids (total vanillic acid 
syringic total acid) with PC2, unlike most hydroxyninnamic acids 
(ferulic acid, caffeic acid and p-coumaric acids) are correlated 
with PCI. This may be because the hydroxybenzoic acids and 
hydroxyninnamic acids belong to different pathway and these 
pathways are not strongly correlated. 

Several authors have shown high and positive correlations 
between total free, conjugated and bound phenolic acids (Verma 
et al., 2008; Fernandez-Orozco et al., 2010). In our study there 
is a strong correlation between total conjugated and united phe- 
nolic. On the contrary, the total free phenolics is less correlated 
with the above (total bound and conjugated phenolics) in con- 
trast to that reported by Fernandez-Orozco et al. (2010). Syringic 
acid is the compound with mayor contribution to the total free 
phenolics, while ferulic acid is the main phenolic compound of 
bound and conjugated phenolics. The low correlation observed 
between the free phenolics with conjugated and bound phe- 
nolics may be due to that the ferulic acid and syringic acid 
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belong to different biosynthetic pathways, as we have described 
above. 

The environmental effect (the growing season and location) 
has been proposed as the main factor that determines phenolic 
acid content. Environmental factors such as solar radiation, tem- 
perature, and irrigation, affect the content of phenolic acids in 
wheat but some phenolics fractions are more susceptible to the 
effect of environment than others (Moore et al., 2006; Fernandez- 
Orozco et al, 2010). Mpofu et al. (2006) has shown that the 
effect of the environment is greater than the effect of genotype 
on hard and spring wheat. Our study corroborates the fact that 
the environment has a strong effect on the content of pheno- 
lic acids, as shows the effect of the factor "Year" in our analysis. 
Conversely, our study shows that although the effect of the envi- 
ronment is high the genotype has a greater effect in many of the 
analyzed phenolics fractions. In fact, the factor "Year" only had 
a greater impact than genotype on the total ferulic acid, total 
p-coumeric, bound p-coumeric and total phenolic. Our study 
shows a greater effect of factor "genotype" that has not been pre- 
viously shown in wheat. This may be due to the higher genetic 
variability of tritordeum in comparison to wheat. This variabil- 
ity gives amphiploid tritrordeum a high value as source of genetic 
variability to improve the antioxidant content of wheat, or as a 
cereal food crop in its own right. 

As shown in Figure 4, the chromosome substitution 
DS2D(2H'^^) increases the content of hydroxyninnamic acids 
(caffeic acid and p-coumaric acids). We also observed an effect on 
the phenolic acid content of the substitution DS1D(1H'^''), since it 
increases the total dimers. This suggest that either the absence of 
chromosome 2H'^'^ or the introgression of chromosome 2D leads 
to increase the content of phenolic compounds. The widely used 
translocation IRS/IBL increases the content of total, conjugated 
and bound ferulic acid, bound p-coumaric acid and total bound 
phenolics in tritordeum. This is the first work where the phenolic 
acid content is associated with the substitution of chromosome 
DSlD(lff'i) and DS2D(2ff*'), and the translocation IRS/lBLin 
tritordeum. Therefore, on chromosome 1 and 2 may be located 
genes involved in hydroxyninnamic acids biosynthesis and/or 
accumulation. 

In conclusion, the high variability of phenolic compounds 
present in tritordeum confers this a high potential for improv- 
ing the quality and quantity of dietary antioxidants in cereals. 
However, in this connection, some limitations are worth noting. 
The transfer of traits from wild relative species to wheat can be 
achieved by uncontrolled translocations which may also intro- 
duce non-desired traits, or by the identification of genes involved 
in the trait and its use for wheat transformation. Thus, although 
tritordeum is very interesting source of genetic variability for phe- 
nolic acid content more in depth understanding of the genes 
that determine the content of phenolic compounds is needed for 
directed improvement of antioxidants in tritordeum or in related 
cereals. 

FUNDING 

This work was supported by the Spanish Comision 
Interministerial de Ciencia y Tecnologia (AGL201 1-22596). 
FP is supported by "Ramon y Cajal" research contract from the 
MINECO (RYC-20 10-07345). 



ACKNOWLEDGMENTS 

The authors thank Carmen Ramirez for her dedication to the 
development and care of the plant material used in this work. 
The authors thank Dr. Paul A. Lazzeri for his critical review of 
the manuscript. 

SUPPLEMENTARY MATERIAL 

The Supplementary Material for this article can be found online 
at: http://www.frontiersin.org/journal/10.3389/fpls.2014.00010/ 
abstract 

Table S1 | ANOVAs tables from the MANOVA analysis of the Model 1 
('variable ~ Line + Year') and IVIodel 2 ('variable ~ Year + IRS/1 BL + 
Substitutions + Cytoplasm + Year x 1RS/1BL + Year x Substitutions -I- 
Year x Cytoplasm + IRS/1 BL x Substitutions + IRS/1 BL x Cytoplasm + 
Substitutions x Cytoplasm'). 

Figure SI | RP-HPLC chromatograms of the pure phenolic compounds 
studied in this work. For this analysis we used pure phenolics from 
SIGMA-ALDRICH. Mix, a mix of the same amount of each phenolic acids. 
FA, ferulic acid; pCA, p-coumaric acid; CA, caffeic acid;VA, vanillic acid; SA, 
syringic acid. 

Figure S2 | Spectra of the pure phenolic acids. 

REFERENCES 

Adorn, K. K., Sorrells, M. E., and Liu, R. H. (2003). Phytochemical profiles and 
antioxidant activity of wheat varieties. /. Agric. Food Chem. 51, 7825-7834. doi: 
10.1021/jf030404l 

Atienza, S. G., Ballesteros, J., Martin, A., and Hornero-Mendez, D. (2007). 
Genetic variability of carotenoid concentration and degree of esterifi- 
cation among tritordeum (x Tritordeum Ascherson et Graebner) and 
durum wheat accessions. /. Agric. Food Chem. 55, 4244-4251. doi: 10.1021/ 
jf070342p 

Beta, T., Nam, S., Dexter, J. E., and Sapirstein, H. D. (2005). Phenolic content and 
antioxidant activity of pearled wheat and roller-milled fractions. Cereal Chem. 
82, 390-393. doi: 10.1094/CC-82-0390 

Castillo, A., Ramirez, M. C, Martin, A. C, Kilian, A., Martin, A., and Atienza, S. 
G. (2013). High-throughput genotyping of wheat-barley amphiploids utilising 
diversity array technology (DArT) . BMC Plant Biology 13:87. doi: 10. 1 1 86/ 147 1 - 
2229-13-87 

de Munter, J. S. L., Hu, F. B., Spiegelman, D., Franz, M., and van Dam, R. M. 
(2007). Whole grain, bran, and germ intake and risk of type 2 diabetes: 
a prospective cohort study and systematic review. PLoS Med. 4:e261. doi: 
1 0. 1 37 l/journal.pmed.004026 1 

Fernandez-Orozco, R., Li, L., Harflett, C., Shewry, R R., and Ward, J. L. (2010). 
Effects of environment and genotype on phenolic acids in wheat in the 
HEALTHGRAIN diversity screen. /. Agric. Food Chem. 58, 9341-9352. doi: 
10.1021/jfl02017s 

Flight, I., and Clifton, P. (2006). Cereal grains and legumes in the prevention of 

coronary heart disease and stroke: a review of the literature. Eur. J. Clin. Nutr. 

60, 1145-1159. doi: 10.1038/sj.ejcn.l602435 
Fox, J., and Weisberg, S. (2011). An R Companion to Applied Regression, 2nd Edn. 

Thousand Oaks, CA: Sage. 
Hothorn, T, Bretz, R, and WestfaU, P. (2008). Simultaneous inference in general 

parametric models. Biom. f. 50, 346-363. doi: 10.1002/bimj.200810425 
Husson, R, Josse, J., Le, S., and Mazet, J. (2013). FactoMineR: Multivariate 

Exploratory Data Analysis and Data Mining with R. R package version 1.25. 

Available online at: http://CRAN.R-project.org/package=FactoMineR 
Irakli, M. N., Samanidou, V. R, Biliaderis, C. G., and Papadoyannis, I. N. (2012). 

Development and validation of an HPLC-method for determination of free and 

bound phenolic acids in cereals after solid-phase extraction. Food Chem. 134, 

1624-1632. doi: 10.1016/j.foodchem.2012.03.046 
Jonnalagadda, S. S., Harnack, L., Liu, R. H., McKeown, N., Seal, C, Liu, S., et 

al. (2011). Putting the whole grain puzzle together: health benefits associated 

with whole grains — summary of American Society for Nutrition 2010 satellite 

symposium./. Nutr. 141, 1011S-1022S. doi: 10. 3945/jn.ll0.132944 
Kaiser, H. F. (1974). An index of factorial simplicity. Psychometrika 39, 31-36. doi: 

10.1007/BF02291575 



www.frontiersin.org 



January 2014 | Volume 5 | Article 10 | 9 



Navas-Lopez et al. 



Phenolic in tritordeum 



Krygier, K., Sosulski, R, and Hogge, L. (1982). Free, esterified, and insoluble-bound 

phenolic acids. 1. extraction and purification procedure. /. Agric. Food Chem. 30, 

330-334. doi: 10.1021/jf00110a028 
Lewin-Koh, N. J., Bivand, R., Pebesma, C. B. E. J., Archer, E., Baddeley, 

A., Bibiko, H.-J., et al. (2012). Maptook: Tools for Reading and Handling 

Spatial Objects. R package version 0.8-14. Available online at: http://CRAN.R- 

project.org/package— maptools 
Li, L., Shewry, R R., and Ward, J. L. (2008). Phenolic acids in wheat varieties in 

the HEALTHGRAIN diversity screen. /. Agric. Food Chem. 56, 9732-9739. doi: 

10.1021/jf801069s 

Marquardt, L., Asp, N. G., Richardson, P., Kamp, J. W., Asp, N. G., Miller Jones, ]., 
et al. (2004). "Whole grain health claims in the United States, United Kingdom 
and Sweden," in Dietary Fibre — Bioactive Carbohydrates in Food and Feed, eds 
J. W. van der Kam, N. G. Asp, J. Miller Jones, and G. Schaafsma (Wageningen: 
Wageningen Academic Publishers), 39-57. 

Martin, A., and Sanchez-Mongelaguna, E. (1982). Cytology and morphology of the 
amphiploid Hordeum chilense x Triticum turgidum conv. Durum. Euphytica 31, 
261-267. doi: 10.1007/BF00028329 

Martin, A., Alvarez, J., Martin, L., Barro, R, and Ballesteros, J. (1999). The develop- 
ment of tritordeum: a novel cereal for food processing. /. Cereal Sri. 30, 85-95. 
doi: 10.1006/icrs.l998.0235 

Martin, A., Martin, L. M., Cabrera, A., Ramirez, M. C, Gimenez, M. J., Rubiales, 
D., et al. ( 1 998) . The potential of Hordeum chilense in breeding triticeae species. 
Triticeae in, 377-386. 

Mattila, P., Pihlava, J.-M., and Hellstrom, J. (2005). Contents of phenolic acids, 
alkyl- and alkenylresorcinols, and avenanthramides in commercial grain prod- 
ucts. /. Agric Food Chem. 53, 8290-8295. doi: 10.1021/jf051437z 

Mellado-Ortega, E., and Hornero-Mendez, D. (2012). Isolation and iden- 
tification of lutein esters, including their regioisomers, in tritordeum 
(x Tritordeum Ascherson et Graebner) grains: evidence for a preferen- 
tial xanthophyll acyltransferase activity. Pood Chem. 135, 1344—1352. doi: 
10.1016/j.foodchem.2012.05.046 

Moore, J., Hao, Z., Zhou, K., Luther, M., Costa, J., and Yu, L. L. (2005). Carotenoid, 
tocopherol, phenolic acid, and antioxidant properties of maryland-grown soft 
wheat. ;. Agric Food Chem. 53, 6649-6657. doi: 10.1021/jf050481b 

Moore, J., Liu, J.-G., Zhou, K., and Yu, L. L. (2006). Effects of genotype and envi- 
ronment on the antioxidant properties of hard winter wheat bran. /. Agric. Food 
Chem. 54, 5313-5322. doi: 10.1021/jf0603811 

Mpofii, A., Sapirstein, H. D., and Beta, T. (2006). Genotype and environmen- 
tal variation in phenolic content, phenolic acid composition, and antioxi- 
dant activity of hard spring wheat. /. Agric. Food Chem. 54, 1265—1270. doi: 
10.1021/jf052683d 



Okarter, N., Liu, C.-S., Sorrells, M. E., and Liu, R, H. (2010). Phytochemical content 

and antioxidant activity of six diverse varieties of whole wheat. Food Chem. 119, 

249-257. doi: 10.1016/i.foodchem.2009.06.021 
R Development Core Team. (2012). R: A Language and Environment for Statistical 

Computing. Vienna, Austria: R Poundation for Statistical Computing. ISBN 3- 

900051-07-0. 

Shewry, P. R., Piironen, V., Lampi, A.-M., Edelmann, M., Kariluoto, S., Nurmi, T., et 
al. (2010). The HEALTHGRAIN wheat diversity screen: effects of genotype and 
environment on phytochemicals and dietary fiber components. /. Agric. Food 
Chem. 58, 9291-9298. doi: 10.1021/jfl 00039b 

Verma, B., Hucl, P., and Chibbar, R. N. (2008). Phenolic content and antioxi- 
dant properties of bran in 51 wheat cultivars. Cereal Chem. 85, 544-549. doi: 
10. 1094/CCHEM-85-4-0544 

Yu, L., and Zhou, K. (2005). Antioxidant properties of bran extracts from 
'Platte' wheat grown at different locations. Food Chem. 90, 311-316. doi: 
10.1016/j.foodchem.2004.04.007 

Zhou, K., Su, L., and Yu, L. L. (2004). Phytochemicals and antioxidant properties 
in wheat bran. /. A^'c. Food Chem. 52,6108-6114. doi: 10.1021/jf049214g 

Zhou, K., Yin, J.-J., and Yu, L. L. (2005). Phenolic acid, tocopherol and carotenoid 
compositions, and antioxidant functions of hard red winter wheat bran. /. 
Agric. Food Chem. 53, 3916-3922. doi: 10.1021/jfD50117c 

Conflict of Interest Statement: The authors declare that the research was con- 
ducted in the absence of any commercial or financial relationships that coiJd be 
construed as a potential conflict of interest. 

Received: 31 October 2013; accepted: 07 January 2014; published online: 30 January 
2014. 

Citation: Navas-Lopez JF, Ostos-Garrido FJ, Castillo A, Martin A, Gimenez MJ 
and Piston F (2014) Phenolic content variability and its chromosome location in 
tritordeum. Front Plant Sci. 5:10. doi: 10.3389/fpls.2014.00010 
This article was submitted to Crop Science and Horticulture, a section of the journal 
Frontiers in Plant Science. 

Copyright © 2014 Navas-Lopez, Ostos-Garrido, Castillo, Martin, Gimenez and 
Piston. This is an open-access article distributed under the terms of the Creative 
Commons Attribution License (CC BY). The use, distribution or reproduction in other 
forums is permitted, provided the original author(s) or licensor are credited and that 
the original publication in this journal is cited, in accordance with accepted academic 
practice. No use, distribution or reproduction is permitted which does not comply with 
these terms. 



Frontiers in Plant Science | Crop Science and Horticulture 



January 2014 | Volume 5 | Article 10 | 10 



